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Abstract 



(^ ■ We have extended recently proposed model of parton distributions in nu- 



cleons to the case of nucleons in nuclei. 
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Recently a simple model for parton distributions in hadrons has been presented 
[Q. They are derived from a spherically symmetric, Gaussian distribution, the width 
of which reflects, via the Heinserberg uncertainty relation, the hadronic size. Two 
distinct parts are distinguished in a hadron: a "bare" hadron (identified with va- 
lence quarks and gluons) and hadronic fluctuations (identified with pions, which are 
later the source of sea partons). The parton density distributions were calculated 
numerically using a Monte Carlo technique, and good agreement with deep inelastic 
scattering data was reported |T[]. 

In this paper we have used the same method to calculate parton distributions 
in nuclei. Preserving the simplicity of the model and using the standard nuclear 
structure we were able to successfully describe the F2{x) in nucleus over the whole 
range of x. The only changes introduced were dictated by the presence of other 
hadrons around the one under investigation in a similar fashion as in our previous 
work on this subject [Q. However, in the present case we were able to describe also 
the region of small x, which shows effects of shadowing and therefore is of particular 
interest to the physics of heavy ion collisions [^ . 

Let us first recollect the main features of the model [|l|, which we shall follow 
as close as possible (in particular the choice of kinematics and notation will be the 
same). Hadron is visualised there as a composition of bare hadron and hadronic fluc- 
tuations. The former are made out of valence quarks and gluons and the latter are 
the source of sea quarks and gluons. They are formed mainly by the pion. One starts 
with the hadron at rest in which frame all partons are supposed to be distributed 
according to a spherically symmetric, Gaussian distributions. Such a form is natural 
because of a large number of interactions binding partons in the hadron. The only 
parameter here is the width of this distribution, the value of which is expected to be 
of the order of a few hundred MeV (both from the Heinserberg uncertainty relation 
applied to the hadron size and from the primordial transverse momentum of partons 
observed in deep inelastic collisions). This parameter encompasses also all pertur- 
bative QCD effects present due to the initial state emission and therefore depends 
on the scale Qq Q The goal of this approach is not so much the full wave function 
of the hadron, as the probability of finding a single parton with the four momentum 
k probed by external current with four-momentum q. Therefore all other partons 
are treated collectively as a single remnant with the four momentum r. Because the 



^Actually, when confronting flavour dependent data the number of parameters is enlarged be- 
cause each quark flavour can have its own width; the same is true for gluons, cf. y. In such 
situation the percentages of each species in the momentum sum rule is also a kind of parameter. 



above prescription provides us only with the three-momentum of the probed parton 
it is assumed that the energy component is equal to the parton (current) mass plus 
a Gaussian variation with the same width as above. It means that the parton can 
be off-shell at the scale Qo and fluctuates with a life-time given by the nucleon radius. 

The reaction takes place in a coordinate system in which the negative z-direction 
points along the current which probes the hadron. One uses the light cone momen- 
tum fraction x of the parton defined as fc+ = xp^ (where p is the four-momentum 
of the hadron). The final four- momentum of the scattered parton denoted by j 
must satisfy the following condition: < j^ < W"^ (where W is invariant mass of 
the hadronic system). This is equivalent to < x < 1. When masses of quarks 
are neglected the same condition must be satisfied by the four-momentum r of the 
remnant: < r^ < W'^. The parton distributions are then calculated by a Monte 
Carlo code 0. The momentum k component of the parton to be probed by current 
with virtuality Qq and four-momentum q is chosen from the Gaussian distribution 
described before. Actually the values of Qg is a free parameter expected to be of 
the order of 1 GeV^. This makes it possible to calculate the four-momenta j and r. 
Events are accepted if the exact kinematical constraints mentioned above are ful- 
filled. In this way one obtains as a result the valence parton distribution fy{x; Ql) 
("bare" hadron) |. The sea parton distribution is given by the pionic component 
of the nucleon: /^(x; Ql) = J ^Uiv, Ql)fpion{x/y; Ql). It is the convolution of the 
pion distribution function in the nucleon, fn{x; Ql), and the parton structure of pion 
fpion{x', Ql). This in turn is obtained from the same Gaussian primodial distribution 
as used for valence partons. The characteristic behaviour of the sea partons is then 
derived from the pion distribution in the nucleon, which was again parametrized 



^In our work we have used a simplified version of model |lj] witli the same distributions for all 
valence quarks and without evolution in Q^. We shall also not address gluon distributions, which 
in Q were assumed to have the same basic Gaussian shape as the valence quarks. The presence 
of gluons will be accounted for only in the momentum sum rule, where part of momentum will 
always be allocated to the gluonic component. 

■^One has to remember, however, that it is, the so called, Bjorken variable x = xbj = 2~, 
which is accessible experimentally whereas |l| starts from the light cone target rest frame variable 
X — xlc — ^ (where p = (Af , 0, 0, 0)) with a fixed resolution Q"^ = Q^. Whenever one is interested 
in some higher Q^, the QCD evolution in the initial state must be performed. In this case the 
X — xbj picked out of the proton at Q^ will be smaller than the xlc — k^ /p'^ one has started 
with at Qq. The remaining energy is radiated into the final state as jets. If the invariant mass 
squared of all jets in the final state is M^, then x^c — ^Bj ( 1 + -q? ) • I^ this case, if one only 

scatters a quark with no further radiation, then x^c — xbj ( 1 + qt ) — ^Bj- This clarifies the 
limits of applicability of the model used here. 



by Gaussian distribution with a smaller width equal to 0.052 GeV reflecting the 
smallness of the pion mass. The total nucleon structure function is then equal to: 
F{x;Ql) = Ux,Ql) + ax;Ql). 



We shall now apply this model to deep inelastic collisions with nuclei, / + A — > 
I' + X. As usual in such cases our aim will be the nuclear structure function F^{xa), 
which shows a characteristic pattern: shadowing at small x, followed by antishad- 
owing ata:;~0.1-^0.3, followed in turn by a very pronounced deep around x ~ 0.7 
and a kind of cumulative effect for x ^ 1. This subject has already long history and 
a vast literature (cf. |^ for review). Our aim is to see if, and under what conditions, 
the model proposed in [^ for free nucleons can be also applied to nuclei. In the 
approach presented here, from many possibilities mentioned in 0], we have selected 
the picture in which collision proceeds on nucleons bounded in nuclei, looking for 
nuclear partonic distributions as a sum of distributions of bound nucleons. The 
corresponding nuclear structure function can be then written as simple convolution 
of nuclear and nucleonic components: 

J- -.— 4 / s . /" , f CLX 

xa 

Here F^ {x) denotes a nucleon structure function inside the nucleus and p^{yA) 
is the nucleon distribution function in the nucleus. Variables xa/A and Va/A are 
the corresponding Bjorken variables for the quarks and nucleons in the nucleus, 
respectively (with nucleonic and nuclear longitudinal momenta given by p'^ and Pa)- 
The presence of nucleus is summarised here by the nucleon distribution function p 
and the simplest approach one could think of would be to use available information 
on it and keep the nucleonic structure functions unchanged. For example, in the 
framework of the relativistic mean fleld theory (RFM) [Q leptons interact with 
nucleons with density p^, which are moving in a constant average scalar Us and 
vector Uv potentials (deflned in the rest frame of the nucleus as Us = —400— MeV 



-F^{xa) = AJdvA l—6{xA~yAx)p\yA)F^{x). (1) 
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and Uv = 300^ MeV, in both cases po = 0.17 fm'^). As a resuh (cf. 0) 
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the Relativistic Fermi Gas approximation the following form of the nuclear density 
function is obtained: 
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The p is nuclear density (which is different from the density of infinite nuclear matter 
Po); Sn = n{p)6 {p^ — [E*{p) + Uv]} denotes nucleon spectral function with n{p) 
being Fermi distribution of nucleon momenta inside the nucleus, p G {0,pf). The 
corresponding Fermi momentum pp = (3/2p7r^)^/^. The nucleon chemical potential 



/i = Ma/ A (which in RMF can be shown to be /x = E*{pf) + Uy) and E*{p) = 
y/p'^ + m*'^ with m* = m — Us- After integration (H) simphfies to 

pM-4 ^1 • (3) 

Here va = Pp/Ep and i/a takes the values given by the inequahty: < {E*p —pp) < 
fii/A < {Ep + pf)- In this way the motion of the nucleon inside the nucleus is 
parametrized here by two parameters: nuclear density p which determines Fermi 
momentum pp and nucleon chemical potential /i. Out of these two, the nucleon 
chemical potential is essentially constant (and equal to /i = 8 MeV), except for a 
few very light nuclei. The nuclear density p due to the finite size of the nucleus 
can vary from p ~ 0.1 fm~^ for light nuclei to p ^ 0.17 fm~'^ for heavy nuclei 
and nuclear matter. Taking p ^ 0.12 fm~^ for the average density of nucleons in 
^^Fe with Pp = 240MeV and using (^) results in dashed curve in Fig. 1, which is 
completely off data. This means that RFM approach to deep inelastic scattering 
on nuclei fails, at least in the form presented here. One could play with nuclear 
parameters, for example by increasing chemical potential p to make the structure of 
deep around x ~ 0.5 -^ 0.7 more pronounced and more similar to the experimental 
data. This would, however, be difficult to justify and the shadowing/antishadowing 
structure seen in data at smaller x would remain unexplained. There are then two 
possibilities: either to keep nucleon structure functions unchanged and try to change 
p by including correlations ^^ and/or cluster structure of nucleons inside the nu- 
cleus |]T^ or to decide to change the nucleonic structure functions. Because the first 



approach will essentially not affect the shadowing (or even antishadowing) region 
of X, we have opted for the second possibility. In what follows we shall therefore 
use, for nuclear structure function, the nuclear density p^ as given by eq. (|^) with 
p = 8 MeV and with Us = Uy = (i.e., m* = m). This is to avoid the possible 
double counting of nuclear interactions, which in our model can be expressed either 
by potentials Uv,s (treated now as free parameters) or by changes of the free nucleon 
structure functions discussed below. Because we have decided (as said above) for 
the latter possibility we have consistently set values of potentials equal to zero. It 
means that we are considering here a model of Fermi gas with modified nucleons. 

Our Monte Carlo calculation for collision on nuclei has been performed in the 
same way as in |^. The only difference was that now the nucleon three-momentum 
is not fixed but chosen from distribution p^{iia) as given by eq. (||), i.e., we account 
for the Fermi motion of nucleons in nucleus. This introduces a kinematical medium 
effect produced by Lorentz transformation of the parton momenta from the nucleonic 



to the nuclear rest frame. This Fermi motion of nucleons also affects invariant mass 
W, four momentum r = p — j and Bjorken variable x = k^/p^ and in this way 
influences the calculated structure function F^{x,p) of the bound nucleon making 
it momentum dependent. Analytically it would mean that one extends our simple 
convolution formula (|l]) replacing F^{x) by momentum dependent F^{x,p) and 
integrates over momentum p, i.e., instead of (|1|) one uses formula 



F2^(^A) = ^ldyAl^,sM 
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The other modifications mentioned before must be made by some suitable changes 
in parameters of the original nucleonic structure functions. Because we are not 
differentiating here between partons of different flavour there are only three such 
parameters, which can be argued to be affected by nuclear medium: the widths of 
the original gaussians (i.e., transverse primordial distributions), the same and equal 
(jq = 0.18 GeV for both valence quarks ("bare" hadron) and for nucleonic pions (sea 
quarks) and the width of the pionic distribution in nucleon equal cr,r = 0.052 GeV. It 
turns out that, in order to get reasonable description of data, it is enough to modify 
only aq (keeping it again the same for both types of quarks) by decreasing it to the 
value (jq = 0.165 GeV and centering pionic energy distribution not at the pion mass 
TTi^ = 0.14 GeV but at rriTr = 0. This can be seen in Fig. 1 where we present results 
of calculations of ratio R{x) = f'z'ix)/ f^{x) for ^^Fe (performed using average value 
of nuclear density p = 0.12 fm~^, pp = 240MeV). The mentioned results are pre- 
sented by open squares. Such a change in the width of initial partonic distribution 
can be naturally explained by the expected swelling of nucleons in nuclear matter 
P, ^ . Smaller spread of momenta caused by the internucleonic interactions corre- 
sponds, due to the uncertainty relation, to bigger spread in coordinate space (the 
"vacuum" in nucleus is obviously different from that outside the nucleus). Although 
the idea is similar to what was proposed a long time ago in p[ the effect of swelling 
is this time modelled not by changing the bound nucleon mass but by changing the 
intrinsic motion of partons inside such nucleon. What concerns the change made in 
the center of the energy distribution of nucleonic pions one can argue, following for 
example |jl3|, that in nuclear matter such pions, interacting strongly with nuclear 



matter and coupled to the nuclear collective modes, seem to have indeed vanishing 
effective mass. In any case such change is essential in getting agreement with data. 
With only swelled nucleon and without decreasing m^ (understood as a position 
of maximum in the energy distribution of nucleonic pions) results in curve denoted 
by open diamonds in Fig. 1. The other parameters were left unchanged (in par- 
ticular, the fraction of momentum carried by the sea quarks is left 7.7% as in [0). 
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Contrary to what has been presented in p[ this time we have also reproduced (at 
least partially) the shadowing/ ant ishadowing region of x. It has to be admited here, 
however, that there is a price in obtaining these results (see open squares in Fig. 
1). The momentum sum rule is now violated by ~ 2%, which in our case means 
that this fraction of nucleonic momenta is taken by gluons (in addition to what they 
already had in free nucleons). If we strictly impose the momentum sum rule (i.e., if 
we assume that there is no additional transfer of the momenta to gluons), it results 
in changing the momentum carried by quarks to 9.7% and spoils the agreement with 
data in the vicinity of the antishadowing region (cf. Fig. 1, open triangles). 

To conclude, we have demonstrated that the simple model of partonic distribu- 
tions in nucleons developed recently in Ref. [^ can also sucessfully describe partonic 
distributions in nuclei by: (z) using standard description of nucleonic Fermi motion, 
{ii) changing by small amount (about 10%) the Gaussian widths of the initial par- 
tonic components and (iii) centering the energy distribution of nucleonic pions on 
777,^ = instead 0.14 GeV. These changes correspond to the valence quark distri- 
butions in nuclei being shifted towards the lower values of x and the sea quark 
distribution being relatively spread out and shifted towards the higher values of x. 
Also nucleons in nucleus allocate a bit more momentum in the gluonic component. 
Altogether, by these simple means a good agreement with data was obtained for the 
whole range of x except of really small values of x < 0.06 where such small changes 
are definitely not sufficient and convolution model is hardly justifiable [§]. We con- 
sider it remarkable that by changing only one parameter describing the primordial 
diffusiveness of partons (in momentum space) in the free nucleon and shifting pi- 
onic distribution in way suggested by nuclear matter calculations we were able to 
describe experimental data in a fairly reasonable way in a very broad region of x. 
Actually, our results are quite robust to changes in initial values of parameters de- 
scribing partonic distributions in free nucleon (for example, a 20% change in the 
quark distribution in pion does not change results presented in Fig. 1). Our anal- 
ysis is obviously simplified, as we have not differentiated between quark fiavours or 
accounted for the presence of gluons only via the momentum sum rule. Aso our 
convolution approach, cf. eq. (|T|) is not suitable for small values of x. But already 
at this level it shows that model [0] is highly suitable for application to nuclear par- 
tonic distributions because of simplicity and cogency of the parametrization chosen 
there. 
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Figure Captions: 



Fig. 1 Results for R{x) = F^{x)/F^{x) for ^^Fe. Data are from p] - full diamonds 
and from ||, [1^ - full circles. See text for explanations. 
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